Nonlinear supratransmission and bistability in the Fermi-Pasta-Ulam model 
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The recently discovered phenomenon of nonlinear supratransmission consists in a sudden increase 
of the amplitude of a transmitted wave triggered by the excitation of nonlinear localized modes of the 
medium. We examine this process for the Fermi-Pasta-Ulam chain, sinusoidally driven at one edge 
and damped at the other. The supratransmission regime occurs for driving frequencies above the 
upper band-edge and originates from direct moving discrete breather creation. We derive approxi- 
mate analytical estimates of the supratransmission threshold, which are in excellent agreement with 
numerics. When analysing the long-time behavior, we discover that, below the supratransmission 
threshold, a conducting stationary state coexists with the insulating one. We explain the bistable 
nature of the energy flux in terms of the excitation of quasi-harmonic extended waves. This leads to 
the analytical calculation of a lower-transmission threshold which is also in reasonable agreement 
with numerical experiments. 

PACS numbers: 05.45.-a; 63.20.Pw; 05.45.Yv 



I. INTRODUCTION 

In a recent series of interesting papers J. Leon and 
coworkers [J [H [J LI discovered that nonlinear chains 
driven at a boundary can propagate energy in the forbid- 
den band gap. Numerical experiments were performed 
for harmonic driving, and the semi-infinite chain ideal- 
ization was simulated by adding damping on the bound- 
ary opposite to driving. In this case, energy transmission 
occurs above a well defined (frequency dependent) crit- 
ical amplitude. This phenomenon has been called non- 
linear supratransmission by the authors, and is charac- 
terized by the propagation of nonlinear localized modes 
(gap solitons) inside the bulk. Several models have been 
considered: sine-Gordon and Klein-Gordon [J, double 
sine-Gordon and Josephson transmission lines 2], Bragg 
media |j| , and an experimental realization has been pro- 
posed for a mechanical system of coupled pendula Q. 
The generic features of the supratransmission instability 
have been described in terms of an evanescent wave desta- 
bilization Q. Moreover, the same process has been de- 
scribed in Ref. [f| for the discrete nonlinear Schrodinger 
equation, suggesting an experimental application to op- 
tical waveguide arrays. 

In this paper we show that the supratransmission phe- 
nomenon ispresent for Fermi-Pasta-Ulam (FPU) nonlin- 
ear chains Q . At variance with all previously considered 
cases, the harmonic driving frequency must lie above the 
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phonon band, since the FPU interparticle potential is 
translationally invariant and, hence, a forbidden lower 
band does not exist (the phonon spectrum begins at zero 
frequency). This entails that the nonlinear modes which 
propagate in the bulk are moving discrete breathers 0- 
Exact static discrete breathers profiles have been pre- 
sented in the literature, but here we use approximate 
analytic expressions for both the low-amplitude solitonic 
case and for the large amplitude situation [8( . This allows 
to perform a study of the instability at the boundary and 
a detailed analysis of the process which leads to the birth 
and the propagation of the discrete breather. By using 
these approximate solutions, we are able to provide ana- 
lytic expressions for the supratransmission critical ampli- 
tudes as a function of the forcing frequencies, which are 
then successfully compared with numerically determined 
values. 

Besides that, we analyse the long-time behavior of the 
system, studying the formation of a stationary state with 
a given energy flux across the chain. The order parame- 
ter of the transition from the insulating to the conduct- 
ing state is, indeed, the average energy flux, which dis- 
plays a jump at the supratrasmission threshold (which 
could then be thought as a sort of non-equilibrium first- 
order transition). We discover that lowering the ampli- 
tude below the threshold, after the stationary state is 
established, does not interrupt trasmission: the conduct- 
ing state survives even at smaller amplitudes and co- 
exists with the insulating state (a sort of bistability is 
present in the system). By further reducing the ampli- 
tude, a threshold appears below which the energy flux 
vanishes without any apparent discontinuity (here we 
have a sort of second-order transition): we develop a 
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theoretical analysis of this new threshold phenomenon, 
which was absent in previous studies. 

The paper is organized as follows. In Section II we 
introduce the model and the equations of motion. Sec- 
tion III deals with the calculation of the energy flux in 
the quasi-linear approximation. Section IV illustrates all 
analytic and numerical results concerning the determi- 
nation of the supratransmission threshold. Section V is 
devoted to the characterization of the stationary states 
and of their bistability. Section VI contains some con- 
clusions. In the Appendix we report, for completeness, a 
calculation of the nonlinear phonon dispersion relation. 



II. THE MODEL 

We consider the Fermi-Pasta- Ulam (FPU) chain [(|, 
which is an extremely well studied nonlinear lattice for 
which a large class of quasi-harmonic and localized solu- 
tions is known. The equations of motion for the so-called 
/3-FPU chain (interparticle potential with a quadratic 
and a quartic term) are 

u n = u n+ i + u n -i - 2u n + (u n +l - u n ) 3 + (lt„_i - u„) 3 , 

where u n stands for the displacement of n-th site in di- 
mensionless units (n = 0, 1, 2 ... , N). All force parame- 
ters have been chosen equal to unity for computational 
convenience. 

To simulate the effect of an impinging wave we impose 
the boundary condition 



u (t) 



A COS LUt 



(2) 



Free boundary conditions are enforced on the other side 
of the chain. 

In order to be able to observe a stationary state in the 
conducting regime we need to steadily remove the energy 
injected in the lattice by the driving force. Thus, we 
damp a certain number of the rightmost sites (typically 
10% of the total) by adding a viscous term —ju n to their 
equations of motion. A convenient indicator to look at is 
the averaged energy flux j = J2 n jn/N, where the local 
flux j n is given by the following formula [tj 

jn = |(«n+ Ml) [u n +l ~ «n + K+l ~ U n ) 3 ] ■ (3) 

Time averages of this quantity are taken in order to char- 
acterize the insulating (zero flux)/conducting (non zero 
flux) state of the system. 



Under the effect of the driving (0), we can look for 
extended quasi-harmonic solutions (nonlinear phonons) 
of the form 



= Acos(kn — cut) 



(4) 



We consider the semi-infinite chain, so that k varies con- 
tinously between and 2ir. The nonlinear dispersion re- 
lation can be found in the rotating wave approximation 
(see e.g. Ref. ^3)- Neglecting higher-order harmonics 
(see the Appendix for details) it reads 



2 A 2 



Lu£(k,A) = 2(1 - cos/c) + 3(1 - coskyA 



(5) 



Thus the nonlinear phonon frequencies range from to 
the upper band-edge u>o(ir, A) > 2. 

If we simply assume that only the resonating phonons 
whose wavenumbers satisfy the condition 



uj = uJo(k, A) 



(6) 



are excited, we can easily estimate the energy flux. Ne- 
glecting, for simplicity, the nonlinear force terms in the 
definition of the flux ©, we have 



J 



-v(k,A)oj 2 A 2 



(7) 



where v is the group velocity as derived from dispersion 
relation (jSJ. This simple result is in very good agree- 
ment with simulations, at least for small enough ampli- 
tudes (see Fig. 0. For A > 0.15 the measured flux is 
larger than the estimate Q, indicating that something 
more complicated occurs in the bulk (possibly, a mul- 
tiphonon transmission) and that higher-order nonlinear 
terms must be taken into account. 
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III. IN-BAND DRIVING: NONLINEAR 
PHONONS 

For illustration, we first discuss the case when the driv- 
ing frequency is located inside the phonon band. Al- 
though trivial, this issue is of importance to better ap- 
preciate the fully nonlinear features described later on. 



FIG. 1: Average energy flux vs. driving amplitude for in- 
band forcing, u) — 1.8, 7 = 5. Data have been averaged over 
10 periods of the driving. The inset is an enlargment of 
the small-amplitude region and the dashed line is the single 
nonlinear phonon approximation @. 
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IV. OUT-BAND DRIVING: 
SUPRATRANSMISSION 

Let us now turn to the more interesting case in which 
the driving frequency lies outside the phonon band, w > 
lvo(tt,0) = 2. In a first series of numerical experiments 
we have initialized the chain at rest and switched on the 
driving at time t = 0. To avoid the formation of sud- 
den shocks ! w e have chosen to increase smoothly the 
amplitude from to the constant value A at a constant 
rate, i.e. 



uq = Acos(cut) 1 — e 



-t/n 



(8) 



where typically we set t\ = 10. 

At variance with the case of in-band forcing, we ob- 
serve a sharp increase of the flux at a given threshold 
amplitude of the driving, see Fig. [3 This phenomenon 
has been denoted as nonlinear supratransmission p] to 
emphasize the role played by nonlinear localized excita- 
tions in triggering the energy flux. 

This situation should be compared with the one of in- 
band driving, shown in Fig.^ where no threshold for con- 
duction exists and the flux increases continuously from 
zero (more or less quadratically in the amplitude). In- 
deed, the main conclusion that can be drawn from the 
previous section is that there cannot be any amplitude 
threshold for energy transmission in the case of in-band 
forcing. Moreover, although at the upper band edge the 
flux vanishes, since it is proportional to the group veloc- 
ity (see formula J7J), it is straightforward to prove that 
it goes to zero with the square root of the distance to 
the band edge frequency. Hence, the sudden jump we 
observe in the out-band case cannot be explained by any 
sort of quasi-linear approximation. 

In the following we investigate the physical origin of 
nonlinear supratransmission, distinguishing the cases of 
small and large amplitudes. 



A. Small amplitudes 

When the driving frequency is only slightly above the 
band (0 < ui — 2 <C 1), one can resort to the contin- 
uum envelope approximation. Since we expect the zone- 
boundary mode k = n to play a major role, we let 



= {-l) n \[^ne iut +r n e 



(9) 



In the rotating wave approximation jjjj and for slowly 
varying tjj n one obtains from the FPU lattice equations 
the nonlinear Schrodinger equation [ip n — > ijj(x,t)) [T^ 

2iu)i> = {uj 2 - A)if) - i> xx - 12V#| 2 , 



(10) 



with the boundary condition ^(0, t) = A. 

The well-known static single-soliton solution of 
Eq. (|10fl corresponds to the family of envelope solitons 
(low-amplitude discrete breathers) 



i„ = <z(— 1)™ cos(wi) sech VG{n — xo)a 



(11) 





Driving amplitude A 



FIG. 2: Average energy flux vs. driving amplitude for out- 
band forcing, uj = 3.5, 7 = 5. Data have been averaged over 
2 10 periods of the driving for a chain of iV = 512 particles. 




lattice site 



FIG. 3: Snapshot of the local energy below the supratrans- 
mission threshold A = 0.15 < A th for u = 2.1 7 = 10. The 
initial condition is an envelope soliton (I1H with xo = —1.8 
(above) and xo = +1.8 (below). 
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with amplitude a = \J {uo 2 — 4) /6. The maximum of the 
soliton shape is fixed by the boundary condition to be 



x = ± 



acosh(a/A) 
ay/6 



(12) 



In this approximation we have two possible solutions: one 
with the maximum outside the chain, which is purely de- 
caying inside the chain (minus sign in l|12|) '). and another 
with the maximum located within the chain (plus sign in 
Overcoming the supratrasmission threshold corre- 
sponds to the disappearence of both solutions. Indeed, 
when the driving amplitude reaches the critical value Ath, 
given by 



lu 2 = 4 + 6A 2 



h i 



(13) 



solution (|llfl ceases to exists. 

We have investigated this issue by simulating the 
lattice dynamics with the initial conditions given by 
Eqs. (fTT|> and l(T2|l . The evolution of the local energy 



ii 2 1 
6 " = ~2 + 2 [V( " Un+1 



,) + V(u n - Un _r)] (14) 



with V(x) = x 2 /2 + x 4 /4, is shown in Fig. 01 The solu- 
tion with the maximum outside the chain (upper figure) 
stabilizes after the emission of a small amount of radi- 
ation (generated by the fact that we have used an ap- 
proximate solution). On the contrary, the other solution 
(lower figure) slowly moves towards the right and, even- 
tually, leaves a localized boundary soliton behind. The 
release of energy to the chain is non stationary and does 
not lead to a conducting state. 

The scenario drastically changes at the supratransmis- 
sion amplitude Ath- The chain starts to conduct: a train 
of travelling envelope solitons is emitted from the left 



boundary (see Fig. Here we should emphasize that 
the envelop soliton solution (f 1 1|) . which is characterized 
by the k = tt carrier wave-number, has a zero group 
velocity. Thus, transmission cannot be realized by such 
envelope solitons. Instead, transmission starts when the 
driving frequency resonates with the frequency of the en- 
velope soliton with carrier wave-number k = tt(N~2) /N, 
next to the 7r-mode. However, as far as we consider a 
large number of oscillators (N = 500), we can still use 
expression (1 1311 for the 7r-mode frequency. 



B. Large amplitudes 

The above soliton solution is valid in the continuum 
envelope limit, and is therefore less and less accurate as 
its amplitude increases. Indeed, if the weakly nonlin- 
ear condition is violated, the width of the envelope soli- 
ton becomes comparable with lattice spacing and, thus, 
one cannot use the continuum envelope approach. Fortu- 
nately, besides the slowly varying envelope soliton solu- 
tion 111|) . an analytic approximate expression exists for 
large amplitude static discrete breather solutions, which 
is obtained from an exact extended plane wave solution 
with "magic" wave-number 2n/3 8] 

u n = a(— 1)" cos [u)B{a)t] cos (J^n ± XqJ , (15) 

if |(7rn/3) ± xo| < tt/2 and u n = otherwise. 
Here xq is defined as follows 

xq = &cos{A/a) , (16) 

where A is the driving amplitude. The breather fre- 




lattice site 



FIG. 4: Snapshot of the local energy at the transmission 
threshold A = 0.253 ^ A th for u = 2.1 7 = 10. The ini- 
tial condition is the envelope soliton [11 1 1 with xq ~ 0. 




FIG. 5: Snapshot of particle displacements u n below the 
supratransmission threshold for a driving frequency ui = 5.12 
and a driving amplitude A — 0.5 < A t h = 2.05. One can 
observe, similarly to the lower Fig. [2] that a moving discrete 
breather appears at the left boundary and propagates inside 
the bulk, leaving behind the static solution. 
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quency <jJb{o) depends on amplitude a as follows 

, ^/3?r 2 (4 + 9a 2 ) 



1.03- 



4AT(s) 



(17) 



where K(s) is the complete elliptic integral of the first 
kind with argument s — 3a/ y/2(9a 2 + 4) and the factor 
1.03 takes into account a reseating of the frequency of the 
"tailed" breather Q (see also |l^|). As previously for 
the case of the envelope soliton solution, we perform a 
numerical experiment where we put initially on the lat- 
tice the breather solution of formula (|15(l . Choosing the 
plus sign in this expression, we do not observe any signifi- 
cant transmission of energy inside the chain. Instead, the 
minus sign causes the appearance of a moving breather, 
which travels inside the chain leaving behind the static 
breather solution with plus sign. Fig. [3] presents this nu- 
merical experiment. 

The static breather solution l|15fl ceases to exist if the 
driving amplitude exceeds the threshold A t h given by the 
resonance condition 



uj = uj B (A th ). 



(18) 



Above this threshold the supratransmission process be- 
gins via the emission of a train of moving breathers 
from the boundary, exactly as it happens in the case of 
small amplitudes. It should be mentioned again that 
the transmission regime is established due to moving dis- 
crete breathers. It has been remarked |8| that discrete 
breathers are characterized by quantized velocities, while 
their frequency is given by the same formula (|17J) . This 
explains why one can use resonance condition i|18|) for the 
static discrete breather solution l|15|l to define the supra- 
transmission threshold in the large amplitude limit. 



C. Supratransmission threshold: numerical test 

To check these predictions, we have performed a nu- 
merical determination of A t h for several values of to, 
starting the chain at rest. This is accomplished by grad- 
ually increasing A and looking for the minimal value A t h 
for which a sizeable energy propagates into the bulk of the 
chain. At early time, the scenario is qualitatively similar 
to the one shown in FigQ] Later on, the interaction of 
nonlinear and quasi-linear modes and their "scattering" 
with the dissipating right boundary establishes a steady 
energy flux into the chain. A conducting steady state, 
which is present also below Ath, will be discussed in Sec- 
tion V in connection with a lower-transmission threshold 

A-th- 

As seen in Fig. [ft] formulae (|18fl [with definition (|17fl ] 
and l|13l) (see the inset) are in excellent agreement with 
simulations for large A > 2 and small A < 1 amplitudes, 
respectively. The accuracy of the analytical estimate in 
formulae (|18f) and l|13|) is of the order of few percents, 
at worst, in the intermediate amplitude range. We don't 




2 4 6 8 10 12 14 16 
Driving Frequency [co] 

FIG. 6: Comparison between analytic estimates and numer- 
ical values of threshold amplitudes vs. the driving frequency. 
Main plot: the full dots are the numerical values of A t h and 
the solid line is a plot of formulae 1)171180 , which are valid for 
large amplitudes. The inset shows an enlargement of the small 
Ath region, in order to illustrate the accuracy of the small- 
amplitude approximation il l -'it (dotted line). The diamonds 
are simulation data for the lower-transmission threshold A~ h 
and the dashed line is formula I22|l . Notice how the latter 
is accurate only for small enough amplitudes (see again the 
inset). 



discuss here the lower curves in Fig. EJ which are related 
to the lower-transmission threshold. 

For comparison, we have checked that the supratrans- 
mission threshold is definitely not associated with the 
quasi-harmonic waves with nonlinear dispersion relation 
(JSJ). If this were the case, the transmission should start 
when the oscillation amplitude reaches the value for 
which the resonance condition to = LUo(k,A) holds. As 
LOo(k, A) is maximal for k = n, we can get the expression 
for the threshold value from the relation to — cjo(ti", A t h), 
i.e. 



lj z = 4- 



(19) 



The amplitude values one obtains from Eq. i|19|) are far 
away from the numerical values and we don't even show 
them in Fig. El This is a further confirmation that supra- 
transmission in the FPU model originates from direct 
discrete breather generation as it happens in the cases 
of discrete sine- Gordon and nonlinear Klein- Gordon lat- 
tices 0. 



V. STATIONARY STATES 

As announced in the Introduction, we have also inves- 
tigated the long-time behavior of the chain. As shown in 
the upper Fig. [7| the time averaged local energy (see for- 
mula dJ) reaches asymptotically a given profile: local 
energy monotonously decreases along the chain as in the 
case of simulations of stationary heat transport with two 
thermal baths |9|. The time-average of the flux © in 
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FIG. 7: Time-averaged local energy (above) and energy flux 
(below) versus lattice position in the case of out-band driving, 
lo = 3.5, A = 1.27. The 50 rightmost particle out of N = 512 
have been damped with 7 = 5. The averages are taken over 
2 10 5 driving periods . 



the stationary state is almost constant along the chain, 
apart from statistical fluctuations and some persistent 
flux oscillations at the left boundary. 

However, as we mentioned above, the value of the sta- 
tionary flux depends on the initial state of the chain. To 
illustrate this effect, let us excite the chain imposing a 
different boundary condition 

u Q = cos(wt) \B(1 - e-*/ Tl ) + (A- B)(l - e^ 2 ) 

(20) 

where T2 3> T\ (in the experiment t-i = 10ri = 100), 
A < A t h and B > Ath- Obviously, both the boundary 
condition (JjJJ and l|20[l lead to the same driving ampli- 
tude A for t ^> r 2 . However, at variance with ©, when 
imposing l|20[l . the istantaneous forcing amplitude over- 
comes the critical amplitude A t h for a time of the or- 
der of T2 , which is enough to establish a stationary flux 
regime. This drastically reduces the transmission thresh- 
old to a value A^ h < A t h, which we denote as lower- 
transmission threshold. This is the first observation of 
this phenomenon, of which we will give a theoretical in- 
terpretation in the following. The numerical determina- 
tion of A^ h versus the driving frequency u> is reported 
with diamonds in Fig. [S] 

In the amplitude interval [^4.7^,^4^], two steady states 
coexist, a conducting state and an insulating one. Each of 
the two steady states can be attained with different initial 
conditions of the chain and different driving pathways. 
For instance, the conducting state is reached when impos- 
ing driving (|20|l . the insulating one when using JSJ). It is 
a typical bistable situation, where two (possibly chaotic) 
attractors coexist in a given control parameter range. 

This behaviour is illustrated in Fig[5] using a different 
simulation method. The average flux is computed after 
changing A stepwise. A back and forth sweep around the 
amplitude interval [A^,^^] reveals the presence of the 



two states. 

A justification of the presence of the lower-trasmission 
threshold can be given in terms of quasi-linear theory. 
This theory leads to dispersion relation © only if one 
restricts to a single right-propagating mode. However, 
due to reflection with the boundary and to mode inter- 
action, both the right-propagating mode and the left- 
propagating one can contribute to the dispersion relation. 
In the Appendix, we derive this more general dispersion 
relation. After introducing the complex mode amplitude 
dfc for the k-th mode, the dispersion relation takes the 
following form 

uj(k) 2 = 2(1 - cos k) + 3(1 - cos kf [\a k \ 2 + 2\a- k \ 2 ] . 

(21) 

In order to fulfill the resonance condition with both the 
right-propagating (a^) and the left-propagating (a_fc) 
mode, their amplitudes must be equal |afc| = |a_fc|. Since 
Lu(k) is maximal for k — ir, the condition for the threshold 
amplitude is 



= 4 + 36 (A 



th) 



(22) 



This analytical estimate (dashed line in Fig. [f)J fits well 
the numerical data only for driving frequencies close to 
the band edge (see the inset). This can be justified by 
taking into account that dispersion relation l|21|l is valid 
only in the weakly nonlinear regime, i.e. mode ampli- 
tudes \ak\ and |a_fc| much smaller than 1. This condi- 
tion certainly applies to the case in which the driving fre- 
quency is close to the band edge, since, then, the thresh- 
old amplitude A7 h is small. When the driving frequency 
is far from the band edge, one has to take into account 
higher-order corrections. The inclusion of the first "satel- 
lite" mode (3fc) produces a lower threshold amplitude, 
but the agreement with numerical data extends only to 
slightly larger amplitudes. To obtain a definitely better 
agreement, one should treat all satellite modes 5fc, 7k, 
etc.. We briefly discuss this aspect in the Appendix. 

From the above considerations, it follows that the 
bistable nature of the energy flux can be explained mak- 
ing reference to the different excitations of the system. 
Indeed, with the system initially at rest, when following 
the driving method J3J), extended quasi-harmonic waves 
cannot be excited. Then, energy flow appears only when 
the driving amplitude reaches the value necessary for lo- 
calized mode excitation. On the other hand, with driving 
(|20|l . the energy flow is initiated by the overcoming of the 
supratransmission threshold and then sustained also by 
extended quasi-harmonic waves. 

It is also possible to give a heuristic argument to ex- 
plain why the transition from non zero to zero flux is 
"continuous" at the lower-transmission threshold A~ h , 
while there is flux jump at the supratransmission thresh- 
old A t h- When the quasi-harmonic waves are already 
excited, reducing the driving amplitude diminishes also 
the number of resonating modes continuously. Hence, the 
flux goes continuously to zero proportionally to this num- 
ber, producing a sort of second-order phase transition, 
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threshold the flux goes to zero continuously. This anal- 
ogy with non-equilibrium phase transitions |17| should 
be further explored. 

Fluctuations in steady states could be analysed to ver- 
ify the possible role played by the Gallavotti-Cohen out- 
of equilibrium fluctuation theorem |l8j . 

The supratransmission phenomenon is quite generic 
and has already been observed experimentally in a chain 
of coupled pendula 0. Also the bistability of conduct- 
ing/insulating states is generic and could be observed 
experimentally in similar conditions. For instance, one 
could apply this theory to micromechanical experiments 
of the type performed by Sievers and coworkers [f9l | . 



Driving amplitude A 



FIG. 8: Energy flux versus driving amplitude for ui — 3, 
N — 512, 7 = 5. The coexistence of the conducting and insu- 
lating regimes is revealed by sweeping the forcing amplitude 
in the range -At/i]. The sweeping direction is indicated 
by the arrows. The analytical values of the two thresholds are 
indicated by the vertical dashed lines. While the prediction 
for the supratransmission threshold Ath is quite good, the one 
for the lower-transmission threshold A~ h overestimates the nu- 
merical value (the prediction becomes better for smaller driv- 
ing frequencies, as shown in Fig. |SJ . In order to show that 
the steady state is already reached for these integration times, 
results for increasing averaging times are displayed with dif- 
ferent symbols. 



when the flux is considered as an order parameter. On 
the contrary, when increasing the driving amplitude with 
the lattice at rest across the supratransmission threshold 
A t h, localized modes are excited, which successively ex- 
cite also extended waves. Hence, a non zero flux is cre- 
ated suddenly from the zero flux state, generating a sort 
of first- order phase transition. 



VI. CONCLUSIONS AND PERSPECTIVES 

We have discussed the supratransmission phenomenon 
for the Fermi-Pasta-Ulam one-dimensional lattice. A 
theory, based on a resonance condition of the driving 
frequency with the typical frequency of localized excita- 
tions (solitons, breathers), gives a good agreement of the 
supratransmission threshold with numerical data. Below 
this threshold two steady states coexist, a conducting and 
an insulating one. For even lower driving amplitudes a 
further transition occurs to a region where only the insu- 
lating state persists: we have called this new phenomenon 
lower-transmission threshold. Imposing a resonance con- 
dition for nonlinear quasi-harmonic waves, we are able to 
derive an analytic expression for the lower-transmission 
threshold amplitude. 

At the supratransmission threshold a jump in the en- 
ergy flux appears. This is reminiscent of a first-order 
phase transition. At variance, at the lower-transition 
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APPENDIX A: NONLINEAR PHONON 
DISPERSION 

In order to derive the nonlinear dispersion relation for 
extended quasi-harmonic waves, let us seek for the solu- 
tions of the equations of motion of the form 



J(uj(p)t+pn) 



+ a + p e~ i{ - u{ - pS>t ~ pnS> 



(AI) 



where u>{p) is the frequency of the p-th mode and a p its 
complex amplitude. Substituting this Fourier expansion 
into the equations of motion, one gets the following infi- 
nite set of algebraic equations for mode amplitudes |16| 

[uj(pf - 2(1 - cosp)] a p = 6 G p qi , q2 a qi a q2 a+ +q2 _ p , 

(A2) 

where 



G 9i,«2 = 7 1 1 + cos (?i + 92) + cos(p - q 2 ) + cos(p - qi ) 



— cosp — cos qi — cos q2 — cos(p — qi — 52)] ■ 

If only a single mode p = k is excited, one gets the 
following dispersion relation 

uo{k) 2 = 2(l-cosfc)+3(I-cosfc) 2 |a fe | 2 , (A3) 

which has been introduced in Eq. JSJ. 
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On the other hand, when both mode k and mode — k 
are excited, one obtains 

u{k) 2 = 2(1 -cos*;) + 3(1- cos fc) 2 [|a fe | 2 + 2|a_ fe | 2 ] , 

(A4) 

which is presented as Eq. (|21[1 in the text. 

As also mentioned in the text, one must sometimes 
consider the excitation of "satellite" modes 3k, 5fc, etc.. 
The inclusion of the 3/c mode produces the addition of 
the following term 

3 [3cos 2 fc- 1 -2cos 3 fc] (|a_ fc | 4 + 2|a_ fc | 2 |a fc | 2 ) , (A5) 



to the r.h.s of Eq. I|A4(1 . This gives the following reso- 
nance condition at k — ir 



^ 2 = 4 + 36(^-J 2 1 + 




where u and A~ h are the driving frequency and lower- 
threshold amplitude, respectively. Since the coefficient 
of the {A^ h ) A term in this relation is always positive, the 
threshold amplitude one obtains is smaller that the one 
derived from Eq. (|22|l in the text. 
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